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The quenched chiral logs are examined on a 16"^ x 28 lattice with Iwasaki gauge action and overlap fermions. 
The pion decay constant f-^ is used to set the lattice spacing, a = 0.200(3) fm. With pion mass as low as 
~ 180 MeV, we see the quenched chiral logs clearly in m^/m and fp, the pseudoscalar decay constant. We 



^ ^ analyze the data to determine how low the pion mass needs to be in order for the quenched one-loop chiral 

pH I perturbation theory (xPT) to apply. With the constrained curve-fitting method, we are able to extract the 

quenched chiral log parameter S together with other low-energy parameters. Only for m„ < 300 MeV do we 

■ obtain a consistent and stable fit with a constant S which we determine to be 0.24(3) (4) (at the chiral scale 

■ 

' ~ ^-^ GeV). By comparing to the 12'^ x 28 lattice, we estimate the finite volume effect to be about 2.7% 

for the smallest pion mass. We also fitted the pion mass to the form for the re-summed cactus diagrams 
and found that its applicable region is extended farther than the range for the one-loop formula, perhaps 
up to ^ 500 — 600 MeV. The scale independent S is determined to be 0.20(3) in this case. We study the 
quenched non-analytic terms in the nucleon mass and find that the coefficient C1/2 in the nucleon mass is 
consistent with the prediction of one- loop xPT. We also obtain the low energy constant L5 from We 
conclude from this study that it is imperative to cover only the range of data with the pion mass less than 
300 MeV in order to examine the chiral behavior of the hadron masses and decay constants in quenched 
QCD and match them with quenched one-loop xPT. 



PACS numbers: 11. 15. Ha, 12.28.Gc, ll.30.Rd 



I. INTRODUCTION 



One of the important goals of lattice QCD is to understand, from first principles, low-energy 
hadron phenomenology as a consequence of chiral symmetry. However, there have been problems 
with regard to chiral symmetry since the advent of the lattice formulation. The lack of chiral 
symmetry in traditional approaches (such as Wilson fermions and staggered fermions) leads to a 
plethora of conceptual and technical difficulties. In particular, due to critical slowing down and 
the existence of exceptional configurations, there is a practical limit on how low in quark (or pion) 
mass one can carry out a Monte Carlo simulation. The existing calculations with these fermions 
are typically limited to pion mass greater than ~ 300 MeV. As such, the extrapolation from the 
current lattice data with pion mass above ~ 300 MeV, sometimes well above, to the physical pion 
mass of 137 MeV is both a theoretical and a practical challenge. The primary concern is whether 
the quark mass region (of the lattice data) is inside the range where xPT applies. If not, one does 
not have a reliable analytic form in quark mass dependence to allow for an unbiased extrapolation. 
Furthermore, most of the quenched calculations are carried out at one-loop order; this limits 
its range of applicability further, to even smaller quark masses. How small does the quark mass 
have to be in order to see the chiral behavior predicted by xPT? Is the strange quark in the radius 
of convergence of xPT? Is there a way to model the chiral extrapolation from well above a 300 MeV 
pion mass without introducing systematic errors? These were some of the issues addressed in the 
panel discussion during the Lattice 2002 conference ^. 

Fortunately, recent progress in lattice fermions has led to formulations (such as domain-wall 
fermions, overlap fermions, and the fixed-point actions) with lattice chiral symmetry at finite cut- 
off. As a consequence, many chiral-symmetry relations P, Q, 0, Q, Ol and the quark propagator 
preserve the same structure as in the continuum. While a lot of progress has been made in 
checking the chiral symmetries in all these three formulations, we find that overlap fermions offer 
several distinct advantages. For example, it has been demonstrated that one can simulate near the 
physical quark- mass region with very gentle critical slowing down P] . Secondly, the overlap- fermion 
inverter can accommodate multiple masses in the calculation of quark propagators 9] and thus 
is particularly well-suited for studying the details of the logarithmically-varying mass dependence 
of physical observables. Thirdly, we find that the O(a^) and 0{'m?a?) errors are small P. Q llo| . 
(There are no 0{a) errors due to chiral symmetry [ll|.) For these reasons, in this paper we use 
overlap fermions in a quenched lattice calculation to study the chiral logs in the pion mass, the 
pseudoscalar decay constant /p, the pion decay constant /vr, and the nucleon mass. We determine 



2 



how small a quark mass or pion mass needs to be for one-loop xPT to be valid. In the present 
work, we only consider hadrons with degenerate quark masses. Quantities involving the strange 
quark such as the kaon mass and decay constant will be studied later. 



II. NUMERICAL DETAILS 

Our calculation is done on a 16^ x 28 lattice with 80 quenched gauge configurations generated 
with the Iwasaki gauge action at /3 = 2.264 and the quark propagators are calculated with 
overlap fermions. The lattice spacing is 0.200(3) fm, as determined from the pion decay constant 
/j^(m7r); the box size is 3.2 fm. 

The massive overlap Dirac op erator [l^ is defined so that at tree-level there is no mass or 
wavefunction renormalization |10l |. 
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D{m)=p+- + {p--)^,e{H). (1) 

Here e{H) = H/ V H'^ is the matrix sign function and Hw is taken to be the hermitian Wilson- 
Dirac operator, i.e. H = Hv/ = 75-Dw- Here is the usual Wilson fermion operator, except 
with a negative mass parameter —p = \/2k — 4 in which < n < 0.25. We take k = 0.19 in our 
calculation which corresponds to p = 1.368. 

Throughout the paper, we shall use lattice units for dimensionful quantities, except that the 
lattice spacing a will be explicit in figures. 



A. Zolotarev Approximation 

The numerical aspects of the calculation have been given elsew^ 



is the adoption of the Zolotarev optimal rational approximation 



To approximate the function , it is suggested \\\ 



nere ,M. The ew ingredient 

n 

UiU^ of the matrix sign function. 



that one should consider the minimization 



of 



|1 - y/xf{x) 



OO 



(2) 



with fix) being approximated by a rational polynomial f{x) G i?™ = pirn — l)/q{m) where 
p[m — 1) and q{m) are polynomials with degree m — 1 and m, respectively. It was proven by 



Zolotarev 



16 1 that G i?^™ is the ||.||oo -optimal rational approximation for the sign 



function on [— 



U [|Amm|, |Amaa;|]- By Way of Zolotarev's theorem [1£, 
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FIG. 1: The residual er{x) of the approximation to the sign function. The soUd Hne is the current 14th 
degree Zolotarev approximation. The dashed and dotted lines are the 14th (IJj and 12th degree Remez 
approximations respectively. 



is an analytic solution 



fix) = A 



where the coefficients can be expressed in terms of Jacobian elliptic functions 

sn^{lK/2m;K) 



ci 



I = 1,2,..., 2m -I, 



(3) 



(4) 



where vT" 



1 - sn'^{lK/2m;K)' 

ly^^l and K is the complete elliptic integral. According to Ref. A is 



determined by the condition 



max [1 — -/xf{x)] = — min [1 — \/xf{2 

c[i, (4^1^)2] c[i,(4^)2] 



(5) 



However, there is a worry when the sign function e{H) is not bounded by unity because of 
slight numerical imprecision. From the Cauchy-Schwartz inequality ||75e(-fr)|| < II75II ||e(-fr)||, the 
overlap operator -D(O) = 1 + j^e^H) with £{H) approximated by e{H) could have superfluous 
unphysical zero modes [2^. In view of this, we fine-tune the parameter A by reducing it a little 
from the value determined from Eq. Q to minimize the positive residual. This doubles the worst 
residuals. 



For the Wilson action kernel Hw = 75-DvF inside the sign function, the largest eigenvalue 
is around 2.43 and is fairly stable over the configuration space. Furthermore, we find that the 
eigenstates of H]y to be projected out become dense at ~ 0.05 — 0.06. Thus, it is sufficient to 
specify \/l — = I ^'""^ I = 80. With lAmaxI = 2.5, this gives lAmml = 0.0315 which is smaller than 
the largest eigenstates (~ 0.05 — 0.06) we project out from the kernel H]y. Since the coefficients of 
the rational polynomial approximation are expressed in terms of the Jacobian elliptic functions, one 
can use as high a degree of the polynomials as the memory can hold in the multi-mass algorithm. In 
practice, we use 14th degree which is a sufficiently good approximation for our calculation. We plot 
the residual €r{x) = -^=—Z\^(x) for x > in Fig.^as a function of x. Here Z\^(x) is the 14th degree 
Zolotarev approximation from Eq. 0. As we can see, in the selected window [0.031, 2.5], the 

approximation is better than 3.3 x lO"^''. This is a good deal better than the earlier attempts which 
approximate g{x) to with the Remez algorithm in which case the resultant approximation to the 
sign function with xg(x^) is not optimal in the range [ lA^T^axl! l-^mml] U [IAtt?,?,^ I i [Am,n..T |]. We plot 
the residual er(x) for the 14th degree jlSj] approximation i?(14) and 12th degree \p\ approximation 
i?(12) with the Remez method in comparison with the 14th degree Zolotarev approximation Z(14). 
We see that the Remez method gives residuals in the order of 10^^ in the relevant range [0.031, 2.5]. 
We note that according to Chebyshev's theorem, the errors of the optimal rational approximation 
f{x) S ij"*'" = p{m)lq(n) have m-|-n-|-2 alternating signs in the range of the approximation. Since 
the Zolotarev approximation to the sign function is optimal, it should satisfy this requirement. In 
other words, there should be 29 alternating maxima and minima (including the endpoints of the 
range) for the approximation with m = 13, n = 14. As seen in Fig. Q this feature is preserved with 
the restriction allowing only positive residuals, as described above. 

After partial fraction expansion, the approximated sign function in rational polynomial xf{x'^) G 

Af=14 



X 

~ x 



r;2 ^ + Qi 

The coefficients q and qi for our 14th degree Zolotarev approximation are listed in Tabled 

In order to improve the approximation of the matrix sign function as well as the convergence 
in the conjugate gradient inversion in the inner do loop P, O, Q], it is desirable to project out 
the lowest several eigenmodes of Hyy. We use the Iwasaki improved gauge action which requires 
that a smaller number of eigenmodes be projected, compared to the Wilson and Liischer-Weisz 
gauge actions [2^. We calculated 140 small eigenvalues A and projected out the lowest 110 - 140 
eigenmodes based on the criterion that the residual of the eigenvalue, defined as || H]y^\ — 



TABLE I: Coefficients for the 14tli degree Zolotarev approximation in Eq. ©. 



1 


0.8371803911435546d- 


-02 


0.4203852824188996d-04 


2 


0.9813674433769438d- 


-02 


0.4230339571009471d-03 


3 


0.1294644813006436d- 


-01 


0.1458829357287322d-02 


4 


0.1831241561050564d- 


-01 


0.3894118671815802d-02 


5 


0.2684553134920763d- 


-01 


0.9481580995276351d-02 


6 


0.4004971283185919d- 


-01 


0.2225210303473774d-01 


7 


0.6031838397347246d- 


-01 


0.5146873615634705d— 01 


8 


0.9155798451341547d- 


-01 


0.1185868564259925d+00 


9 


0.1406107013842408d- 


hOO 


0.2742893835909246d+00 


10 


0.2211980404641135d+00 


0.6437234073136943d+00 


11 


0.3673892837229880d- 


hOO 


0.1567367648339766d+01 


12 


0.6933239005020095d- 


hOO 


0.4183844803034055d+01 


13 


0.1812368256185377d- 


hOl 


0.1442795672202632d+02 


14 


0.1555827970041009d- 


h02 


0.1451886134043600d+03 



A^'a II / II X\/dim.{Hy[/) ||, be less then 10^''. The next smahest eigenvalue is ^ 0.05 - 0.06, which 
is well within the window where the approximation for the sign function is better than 3.3 x 10^^''. 
The residuals of both the inner and the outer do loops for the inversion of the massive overlap 
Dirac operator in Eq. are at the level of 10"''. 

III. PION MASS AND ZERO MODES 

We look at the pion mass as calculated from the pseudoscalar correlator {J2gP{^j (0)) where 
P = — D/2)tp. The pion masses m^r obtained from this correlator are listed in Table ITTl and 

plotted in Fig. [21 as a function of the quark mass m. Our smallest pion mass is 182(8) MeV. 

Consider the contributions of the zero modes to the meson correlator of the interpolation field 



M = '0r(i - L>/2)V', 




+2 E 



i=0,A>0 



(tr(V^l(x)75r^.(x))tr(^|(0)f75V.A(0))) - 

^2 _j_ 
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FIG. 2: The pion mass as calculated from the {P P) correlator as a function of the bare quark mass m from 
0.014 to 0.8. The insert is for the results with small quark masses up to m — 0.1. 

where T is the spinorial matrix and f = ztF^^. For the pseudoscalar case with F = 75, the 
zero modes contribute to both the direct and cross terms in the correlator. They contribute 
also to the cross terms in the ■A^4:{x.,t)P{d)) and Ai{x,t)A4^{Q)) correlators, where = 
^/'i7475(l — D/2)ij), but not to the direct terms. In view of the worry that the pseudoscalar 



Id, 



particularly at small volumes, we 



correlator may be contaminated by the zero modes 
look at A4^{x^t)Ai{Q)) to check if the masses calculated are the same as those from the {PP) 
correlator. We plot the ratios of the pion masses from the {P P) correlator and the (^4^44) 
correlator in Fig. El We see that the ratios are basically unity within errors all the way to the 
smallest quark mass. This suggests that at this volume, i.e. V = (3.2 fm)^, there is no detectable 
contamination due to the zero modes. However, this is not so at smaller volumes. We will show 
results on the 12^ x 28 lattice and discuss the volume dependence in Sec. 



It has been shown that the correlator {PP) — {SS) can get rid of the zero mode contribution [22 1 
for both the direct and cross terms. This is a good procedure for medium quark masses. For small 
quark masses, there is a complication due to the quenched ry'-vr ghost state contribution in the 
(SS) 23] which has a mass close to two times the pion mass. We have fitted the {PP) — {SS) 
correlator with the inclusion of the ghost state which has a negative weight. We show in Fig. 0] 
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TABLE II: Pion mass obtained from the {P P) correlator on a 16"^ x 28 lattice. The pion decay constant 
is used to set the lattice spacing, a = 0.200(3) fm. 



m 




m„ (MeV) 




0.80000 


1.3106(42) 


1293(20) 


0.68 


0.60000 


1.1069(38) 


1092(17) 


0.60 


0.40000 


0.8740(30) 


862 (13) 


0.56 


0.32200 


0.7743(22) 


764 (12) 


0.54 


0.26833 


0.7017(24) 


692 (11) 


0.57 


0.22633 


0.6416(25) 


633 (10) 


0.63 


0.18783 


0.5829(28) 


575 (9) 


0.72 


0.15633 


0.5312(31) 


524 (8) 


0.77 


0.12950 


0.4840(32) 


478 (8) 


0.75 


0.10850 


0.4438(36) 


438 (7) 


0.70 


0.08983 


0.4050(38) 


400 (7) 


0.64 


0.07583 


0.3748(38) 


370 (7) 


0.62 


0.06417 


0.3470(40) 


342 (6) 


0.58 


0.05367 


0.3200(43) 


316 (6) 


0.53 


0.04433 


0.2940(42) 


290 (6) 


0.55 


0.03617 


0.2693(48) 


266 (6) 


0.52 


0.03033 


0.2501(54) 


247 (6) 


0.56 


0.02567 


0.2332(56) 


230 (7) 


0.48 


0.02333 


0.2244(59) 


221 (7) 


0.58 


0.02100 


0.2151(61) 


212 (7) 


0.62 


0.01867 


0.2055(61) 


203 (7) 


0.60 


0.01750 


0.2005(65) 


198 (7) 


0.64 


0.01633 


0.1953(69) 


193 (7) 


0.68 


0.01400 


0.1844(71) 


182 (8) 


0.72 



the ratio of pion masses determined from the {P P) and the (P P) — {S S) correlators. We note 
that the errors from the {P P) — {S S) are larger for m,r < 0.4438(36) due to the presence of 
the ghost state which is close to the ground state pion mass. It is not clear to what extend the 
additional complication due to the ghost state is disentangled from that of the zero modes. In any 
case, we find that the pion masses obtained are consistent with those from the {P P) and {A4 A4) 
correlators. Their central values differ by 2.5% at most. We will come back to this issue in the 
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FIG. 3: Ratio of pion masses as calculated from the {P P) and {A4 A4) correlators as a function of the 
obtained from the {A4 A4) propagator. 



next section when we fit the quenched chiral log in m^. 



IV. CHIRAL LOG IN to. 
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It is predicted that in the quenched approximation, there are quenched chiral logs 
arising from the hairpin diagrams; the flavor-singlet pseudoscalar correlator does not yield a full 
ij' mass due to the absence of dynamical fermion loops. Instead, the would-be r]' propagator gives 
a double pole of the Goldstone boson. The predicted form for the one-loop formula in is the 



following 
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Am{l - 5[ln(^m/A^) + 1]} + A^m'^ll + 2ln{Am/Al)] + Bm^ 



(8) 



where the coefficients A, A^ and B are given in terms of the parameters S, /, Q!$ and 2a% — Q5 in 
the quenched effective chiral lagrangian, i.e. 



A 



Ar 



B 



3(4^' 
(2^8 - 05)^2 
(4vr/)2 ■ 



(9) 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 

ma 

FIG. 4: Ratio of pion masses as calculated from the {P P) and {P P) — {S S) correlators as a function of 
the ma. 

and is an unphysical cutoff scale which is usually taken in xPT to be in the range of 0.5- 
1.0 GeV. The parameter a$ is the singlet coupling in the quenched theory. The quenched chiral 
log parameter 6 is 



(10) 



where mg ~ 870 MeV from the Veneziano model of the r]' mass. From this, one estimates that 
6 = 0.183. 

There are several calculations to extract the chiral log parameter 5 122, 28, 2£, 3C, 31, 32, 
Q, Q, Q. Save for Refs. 0, Q, most of the calculations obtained values for 6 in the range 
0.06 — 0.1 which are smaller than the expectation. We have summarized these results in Table HiH 
The present simulation adopts the overlap fermion action which has exact chiral symmetry on the 
lattice, and a highly accurate numerical approach (the matrix sign function is approximated to 
better than 10^^). This should be a valuable tool to probe the relevant range of applicability of 
the one-loop with pion mass as low as ~ 180 MeV on a fairly large volume (3.2 fm). We 

plot in Fig. El the {PP) results for m^/m as a function of the bare quark mass m. We see clearly 
that below m = 0.1295 (m^r ~ 480 MeV), there is a dramatic divergent behavior. The ratio of the 
highest value at m = 0.014 (m^ ~ 180 MeV) to that at the trough at m = 0.1295 (m^ ~ 480 MeV) 
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is ~ 1.33. Thus, the percentage excess, of the highest value at the minimum pion mass m-7r,mm 
relative to that at the trough at mj^^tTou is 33%. Most of the calculations studying the quenched 
chiral log stopped above ~ 300 MeV; thus, they seem to be seeing the onset of this behavior 
with corresponding excess ~ 6% — 11%. In Table HTll we summarize the attempts made by various 
groups to calculate this quenched chiral log from m^/m with different actions and lattices. Except 
for the present work, all the results are fitted to Eq. (jH)) without the term. 

We see from Table IIIII that the results are not consistent with each other. There is a large 
spread in 5 and mj^^trou which invokes many issues to ponder: To what extent does the fermion 
action matter? Is xPT valid in the range of pion masses calculated? If so, is the one-loop formula 
sufficient and how does one judge it? The parameters A, 6, A-^, and B are not linearly independent 
in Eq. Q; how does one fit them? Practically all the calculations so far ignored A^', does it change 
the results much if included? How large are the finite volume effects for the range of pion masses 
in the calculation? (As will be demonstrated in Section the pion masses calculated from the 
{PP) correlator in small volumes are contaminated by the zero modes even up to fairly large quark 
masses, i.e. the strange quark region. This could, ctffect the position of ff^jT: ^trou 

and the result for 

5.) How does one control the finite volume and zero mode effects? 

Before addressing these questions, we first point out that there is a fundamental problem in the 
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FIG. 5: The ratio m^/m as a function of the quark mass m. 
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TABLE III: Summary of the quenched chiral log calculations in with different fermion actions and ranges 
of quark masses. niTr^min is the minimum pion mass used in the fitting. m,r,trou is the pion mass at the 
trough of rn^/m. £ is the percentage of excess of the highest m^/m at mTr,min relative to the corresponding 
trough at niT.^trou- 



Group 


Fermion 


L a 




^TT.trou 


£ 


S 






(fm) 


(MeV) 


(MeV) 






CP-PACS [28] 


Wilson 


3.2 


300 


- 550 


^ 8% 


0.060(10) - 0.106(5) 


Fermilab [29] 


Wilson- Glover 


2.7 


290 


- 400 


^ 7% 


0.073(20) 


QCDSF [30] 


Wilson- Glover 


1.5-2.3 


400 


- 800 


^ 11% 


0.14(2) 


MILC [31] 


Staggered 


2.6 


320 


- 700 


^ 8% 


0.061(3)) 


Kim k Ohta [32] 


Staggered 


2.6 


210 


- 570 


^ 8% 


0.057(13) 


Blum et al. [22] 


Domain Wall 


3.2 


390 


- 740 


^ 11% 


0.07(4) 


RGB [33] 


Domain Wall 


2.4 


285 


- 660 


~ 6% 


0.107(38) 


BGR [24j 


fixed point 


2.6 


210 


- 870 r 


- 15-20% 


0.17(2) 


BGR [34] 


Ghirally Improved 


2.4 


240 


~ 930 r 


- 15-20% 


0.18(3) 


Ghiu & Hsieh [35] 


Overlap 


1.5 


440 


- 770 


~ 11% 


- 0.2 


Present work 


Overlap 


3.2 


182 


~ 480 


~ 33% 


0.236(33) (37) (one loop) 



0.20(3) (cactus) 



way the lattice data are usually fitted to the form predicted by the xPT. The common practice is 
to fit the lattice data to the result from a high pion mass down to cover the range of calculated 
masses and then extrapolate to the physical pion mass. This can force a fit but cannot answer the 
question as to whether the one- loop chiral formula is applicable for the range of masses considered. 
It is possible that the range of masses will require higher orders in xPT, or even that the mass 
range is outside the radius of convergence of xPT. Then one does not know the precise form to 
fit. Extrapolating down from the pion mass as high as 1 GeV is simply diametrically opposed to 
the philosophy of which is based on the expansion of small quark masses and momenta. 
So, how does one proceed? To be consistent with the yPT approach, we start from the smallest 



36l | . Since we have come down to the pion 



quark mass and ask where the one-loop formula applies 
mass as low as ~ 180 MeV, we will assume that it is in the range where one-loop xPT is valid. We 
fit the one-loop formula with the five smallest pion masses and then incrementally add higher pion 
masses one by one. If the formula is valid for the range of masses, then the parameters will remain 
constant, i.e. they are truly the low energy constants that we seek to evaluate. On the other hand, 
if the parameters begin to change as one includes more higher masses, we take it as an indication 
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that the one-loop formula starts to lose its applicability. This could be due to the fact that higher 
orders are needed or that the mass is simply outside the range of validity of xPT. 



TABLE IV: The low-energy chiral parameters Ci and Cil (and C2 and C2l) are fitted from Eq. with 
a maximum of 17 quark masses in correlated and constrained fits for several ranges of quark masses with 
the mininmm quark mass at 0.0140 (which corresponds to the smallest pion mass at 182(8) MeV). These 
arc then used to determine A, S, Aa, and B from Eq. (|12|) . The parameters A, 6, and B vary with the chiral 
scale, and are evaluated at A^^ — 0.8. 



^max 




Ci 




A 


S 


Acy 


B 


xVdof 


0.02100 


0.2151(61) 


0.82(17) 


-0.368(64) 


1.51(27) 


0.244(40) 


- 


- 


0.48 


0.02333 


0.2244(59) 


0.84(16) 


-0.362(52) 


1.51(24) 


0.239(36) 


- 


- 


0.49 


0.02567 


0.2332(56) 


0.83(15) 


-0.353(51) 


1.48(23) 


0.238(35) 


- 


- 


0.46 


0.03033 


0.2501(54) 


0.84(14) 


-0.350(56) 


1.48(23) 


0.236(35) 


- 


- 


0.44 


0.03617 


0.2693(48) 


0.85(15) 


-0.349(50) 


1.49(22) 
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We mentioned earlier that in the fitting of the one-loop formula, the parameters A, 5, Ay., and 
B are not linearly independent. We can rewrite Eq. Q as 

= Cim + Ci Lfn ln(m) -|- C2m^ -|- C2 ln(m) (11) 

where 

Ci = A[l-S{lnA-2lnAy + l)], 
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FIG. 6: Fitted A, Aa, and B as a, function of mT^_max, the maximum pion mass of the fit, for — 0.8. The 
minimum pion mass is at 0.1844(71) (182(8) MeV). 

CiL = -A6, 

C2 = Aa{l + 2lnA- AlnA^) + B, 
C2L = 2A^, (12) 

are the independent parameters. Therefore, we have one redundant variable which we take to be 



351 with 



the cutoff scale A^. We adopt a Bayesian-based constrained curve- fitting method 
an adaptive procedure to obtain priors sequentially according to the relative importance of the 
parameters. The details of the constrained curve- fitting algorithm for the hadron masses will be 



given elsewhere 



39|. In each of the quark mass ranges, the priors for Ci, Cil, C2, and C2L 



are obtained from the corresponding {A4^A4) correlators. After the C parameters are determined 
from the fit, we solve Eq. H12|) to obtain A, 5, A^, and B for a given A^. We list their results 
evaluated at A-^ = 0.8 together with Ci and Cn in Table HVI We also plot A, A^, B, and 6 as 
a function of the maximum pion mass m-j^^max of the fitting range in Figs. IHl and [Tj We should 
mention that due to the limited number of gauge configurations, we are not able to carry out 
correlated fits with a reasonable for more than 17 quark masses. Therefore, for longer ranges, 
we choose a maximum of 17 representative masses to cover the range of the fit. For nrimax less than 
0.08983 (m^ = 0.4050(38) or 400(7) MeV), we find that the data are dominated by di and Ci. 
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The contributions from C2 and C2 l cancel each other and the combined contribution is neghgible 
compared to those of Cn and Ci. This is illustrated in Fig. El for a typical fit. As a result, it is 
practically impossible to obtain a reliable fit for C2 and C2L. We decided to fit these low mass data 
using a prior for C2L which is limited to less than 2 during the fitting procedure for the (744^44) 
correlators. We tried using the priors without such a limit and also tried dropping C2 and C2 l all 
together in these low mass ranges and found that the fitted Ci and Cn are not changed within 
errors. The central values for 6 are changed only by 1 — 2%. In view of this inability to fit C2 
and C2L, we do not c[uot6 tli6 fitted. cLiid B below TfiYnax — 

0.08983 (m^ = 0.4050(38)) where 
the errors are larger than their central values. Hopefully, with better statistics, one will be able to 
extract the parameters Aa and B in the range m-,^ < 400(7) MeV. 

We see in Fig. El that values for A are fairly stable throughout the range of the fit. A^ and B 
are also stable in the range where rrimax ^ 0.08983 and A^ is consistent with zero in this range. On 
the contrary, 5 in Fig. [Tlis constant up to rrimax ~ 0.044 which corresponds to m,r,maa: ~ 300 MeV. 
Beyond that, it has a sharp drop. We take the fact that 5 is sensitive to the range of masses as an 

0.35 I < < < < < < < < 1 




FIG. 7: The filled circles are the quenched chiral 6 plotted as a function of ruj^^rnax, the maximum pion 
mass of the fitting range, with the minimmn pion mass fixed at 0.1844(71) (182(8) MeV). The open circles 
(filled squares) are those fitted downward from a maximum pion mass fixed at 0.8740(30) (862(13) MeV) 
(0.6414(21) (633(10) MeV)) and plotted as a function of mTr.min, the minimum pion mass of the fitting range. 
These results are at = 0.8. 
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FIG. 8: Fractional contributions for different terms in Eq. The input of the C parameters are from 

the fit with nimax = 0.03617. Note the combined contribution Cim^ +CiLm^ In(TO) denoted by the dashed 
line is very small. 

indication that the one-loop chiral formula in Eq. ((HJ is valid only up to m-,^ ~ 300 MeV. Below 
this value, the parameters in the formula do not exhibit mass dependence and thus they are the 
low energy constants of the effective chiral theory at a certain chiral scale (e.g. A^^ = 0.8). 

Next, we simulate what has been commonly done in lattice calculations. We fix the maximum 
pion mass at = 0.8740(30) (862(13) MeV) and incrementally include lower pion masses to fit 
Eq. (|1H) . We find that no matter how many low pion masses are included, the results are not 
compatible with those starting from rri-j^^rnin = 0.1844(71) (182(8) MeV) and ending at TTIti max ^ 
300 MeV. We illustrate this by plotting, in Fig. [71 the fitted 5 (open circles) as a function of 
the minimum pion mass m-n-^min for the range of the fit. We see that when pion masses down to 
300—400 MeV are included, the fitted 5 is in the range of ~ 0.15 which is consistent with those found 
in Table HTn with a similar mass range and can explain why most of the calculations end up with a 



small 5. The only exception is Ref. 



35l | which obtains a large 6 ~ 0.2 with 771,7^^. 



440 MeV. We 



speculate that this may be due to the small volume used {L = 1.5 fm). As will be explored further 
in Section we find that the zero modes make a sizable contribution to pion masses as high as 
800 MeV on a lattice with L = 2.4 fm, and likely higher at smaller volumes. It would be sensible 
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to check if this is the case by comparing pion masses obtained from the {PP) correlator to those 
from (744^44) as will be done in Section [VJ 

In view of the fact that the rise of /m left of the trough at 171^^ ~ 480 MeV is mostly due to 
the chiral log while the rise right of the trough is mostly due to the term, we first fit the points 
around the trough (i.e. from m = 0.07583 to 0.22633) with the form in Eq. but without the 
Aa term and then extend the range by including smaller quark masses. The results of A, 5 and B 
turn out to be very close to the ones starting from the maximum pion mass at ttItt = 0.8740(30) 
(862(13) MeV). We show the results of 6 (labeled 5trou) as a function of the minimum pion mass 
m-K,niin in Fig. [3 and find them almost the same as 6min at small quark masses and similarly smaller 
than those where the pion masses are restricted to lower than ~ 300 MeV. 
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FIG. 9: The fit to TO^/m as a function of the quark mass m. The solid line is the forward fitting ending at 
= 0.269(5) (to = 0.03617). The dash Une is the backward fitting from = 0.641(2) (m = 0.22633) to 
the smallest quark mass. The insert is a blowup for the small quark mass region. 



For comparison, we show in Fig. |^ the fit of m^/m for the forward fitting ending at = 
0.269(5) (m = 0.03617) (sohd line) and the backward fitting from = 0.641(2) (m = 0.22633) 
to the smallest quark mass (dash line). 

As one can see from Eq. ()12() . the parameters A, 6, and B depend logarithmically on A^. In 
Fig. ^1 we show how A and 6 vary as a function of A^. The Ci and Cil are taken from the fit for 
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rrimax = 0.03617 (m^ = 0.2693(48) or 266(6) MeV) in Table |TV| Incidentally, for a given Ci and 
CiL, there is a maximal value 

(for which 6 = 1 and ^ = Ci) beyond which no solution exists. But this is well beyond the range of 
applicability of Eq. 0. We should point out that the fact A runs with In is a peculiar feature of 
quenched QCD. This is not the case in full QCD where A is basically the quark condensate where 





0.6 0.65 0.7 0.75 0.8 0.85 0.9 0.95 1 



FIG. 10: The running of 6 and A and as a function of according to the one-loop formula. The input of 
Ci and CiL are for the fit with rrimax = 0.03617. 
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there is no quenched chiral S Q]. From Eq. we see that the combination A5 is independent 
of the chiral scale. This makes S also run with In A^. This is not physical, but rather is due to the 
fact that Eq. Q contains only the leading term in the expansion of 6 j4Q|. To see this, consider 
the example where the leading logarithm is re-summed for the "cactus diagrams" j^| (one is not 
concerned with Aa and B here) which leads to 

ml = m^ (AA"^)^ . (14) 

We see that the m and A-^ dependence of S are separable which implies that 5 does not run with 
A^, and the unphysical A^^, can be absorbed in A with 



which gives 



A cx A- (15) 



^ = (16) 

d In A2 ^ ^ 

A. 



This is the same as obtained from Eq. Q when one requires that be independent of A^ and 
assumes that 5 does not depend on A-^. 

We shall report our final results based on the fit of pion mass from 0.184(7) (m = 0.014) to 
0.269(5) (m = 0.03617) with 9 quark masses. This choice is based on the observation that 5 begins 
to decrease beyond rrimax = 0.03617 {mT^max = 0.269(5)). Although this happens to coincide with 
the minimum in x^/dof in Table HVl we should stress that a minimal x^/dof is not and should 
not be the criterion for such a choice. As we pointed out earlier, the primary difficulty in chiral 
extrapolation is the uncertainty about the valid range of one-loop formula in xPT. A /dof less 
than unity does not imply that the fitting formula is correct. One way to judge is the criterion 
that we adopt, namely the one- loop formula is not valid when the fitted parameters do not remain 
constant beyond certain rriT^max- 

As we said earlier, we are not able to obtain a reliable fit for A^, and B at these low masses. As 
for the finite volume systematic errors, we note that the finite volume correction for the lowest 
is estimated to be 5.4% and the second lowest is estimated to be 3.6%. (This will be explained 
in more detail in Section O) We thus estimate the systematic errors due to the finite volume to be 
at the few percent level which are smaller than our current statistical errors. We also dropped the 
last quark mass and fitted the 8 quark masses from m = 0.01633 (rriTj = 0.195(7)) to m = 0.03617 
(m^r = 0.269(5)) and found that A and 5 change only about 1%, which are much smaller than the 
statistical errors. 
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We quote our final results for the low-energy parameters of the quenched one-loop chiral per- 
turbation theory at = 0.8 GeV close to the p mass in Table [V] In view of the unphysical 
dependence of 6 on at the the order specified in Eq. Q, we shall estimate the systematic error 
to be the average deviation from the value at 0.8 GeV of the values at A-^^ = 1 GeV and 0.6 GeV. 
We see that the central value of 5 is somewhat larger than the phenomenological value of 0.183 in 
Eq. (|10|) . but they are consistent within errors. 

TABLE V: The low energy parameters in the quenched chiral lagrangian are given as a fit of Eq. from 
m = 0.014 (m^ = 182(8) MeV) to m = 0.03617 (m^ = 266(6) MeV) at = 0.8 GeV. The second error in 
5 is the estimated systematic error determined as the average deviation at A^ = 1 and 0.6 GeV from the 
value at 0.8 GeV. 
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FIG. 11: The filled circles are the scale independent chiral 5 plotted as a function of mT^^max, the maximum 
pion mass of the fitting range, with the minimum pion mass fixed at 0.1844(71) (182(8) MeV). The same 
are for the triangles {63-) except the starting minimum pion mass is at 0.2055(61) (203(7) MeV. The filled 
squares are those fitted downward from a maximum pion mass fixed at 0.6416(25) (633(10) MeV) and plotted 
as a function of m7r,mirt, the minimum pion mass of the fitting range. 
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FIG. 12: The circles/triangles are fitted A as a function of mTr,„iax starting from m — 0.014/0.01867. The 
squares are those fitted downward starting from m = 0.22633 as a function of mT^^rnin- 




A. Fitting of Re-summed Cactus Diagrams 

Since at small quark masses, the higher order results wi th ( dlnjm))^ corrections are important, 
it is pertinent to fit the re-summed cactus diagram form |24| for the chiral log which leads to a 
scale independent 6. Thus, we shall fit with the power form 



Am^+^ + Bm . 



(17) 



Similar to fitting the one- loop formula in Eq. pi() . we start from the smallest quark mass 
and fit upward first. Since at very small quark masses the B contribution is very small, we put 
a very weak constraint at 1.8 it 1 which covers the range of B as fitted in one-loop formula in 
Table IIVI The resultant 6, A and B are presented in Figs. ^2 ^1 and El We see that, in 
contrast to the sharp transition at ra-j^^raax ~ 300 MeV in the one-loop formula fit, the transition 
for 5 is smoother. The behavior for A and B are similarly changed. This is to be expected, 
since the power form takes into account the re-summed cactus diagrams to all orders. As such, 
it should extend the range of applicability of the quenched xPT to beyond ~ 300 MeV for 
the one-loop case, provided that it dominates the higher loop corrections. From the fitted results, 
it appears that the power form gives a stable fit to the trough of m^/m and is perhaps valid 
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FIG. 13: The circles/triangles are fitted S as a function of rn.^^„iax starting from m — 0.014/0.01867). The 
squares are those fitted downward starting from m — 0.22633 as a function of fn-jj^rnin- 

to ttItt ~ 500 — 600 MeV. Above m-,^ ~ 640 MeV, all the parameters {6, A and B) are different 
from those below 

iT^TT ^max ^ 300 — 400 MeV by two sigmas. We shall Quote the result of the scale 
independent 6 at m = 0.08083(m^ = 400(7) MeV) where 5 = 0.20(3). We have fitted also the pion 
mass from the {P P — S S) correlator and find that the changes of 5 and A are small, at a few 
percent level. In the range between m,r = 0.2501(54) and 0.4840(32), the change is at the 1% level. 
For example, 6 = 0.20(4) at m = 0.08083(m^ = 400(7) MeV). 

We also fitted the points around the trough (m = 0.07583 to 0.22633) and then incrementally 
added smaller quark masses in the fit, as we did for the one-loop formula. The results are plotted 
in Figs. ^2^1 and^las a function of m^^^in and labeled 5min, Amin and Bmin- We see again that 
they do not reproduce those parameters which are obtained from fitting the range from m = 0.014 
to 0.08983 (m^ = 0.4050(38)). 

As we are concerned about the finite volume effect (this will be addressed in Sec. ^ for small 
quark masses, we drop the three smallest quark masses and then fit the pion mass upward. We plot 
the results in Figs. 1111 1121 and as a function of 771,^, max and labeled 5^-,A^^ and i?3-. We see 
that in most cases they are only a few percent different from those {5max^ Amax and Bmax) without 
dropping the last three quark masses. We have also performed the same test for the one-loop log 
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FIG. 14: The filled circles are the chiral S {6max) plotted as a function of 

TTi-ji^max 1 the maximum pion mass 

of the fitting range, with the minimum pion mass fixed at 0.1844(71) (182(8) MeV). The same are for the 
open circles {63^) except the starting minimum pion mass is at 0.2055(61) (203(7) MeV. 

formula fit and find that the parameters are also not changed much. We plot in Fig. E]the results 
of (^3- from this fit together with 6max as obatined from the one-loop log fit without dropping the 
last 3 points. We see that the characteristic transtion around m^r ~ 300 MeV is still visible. 



V. FINITE VOLUME EFFECTS 



Intertwined in the discussion of the validity of the one-loop xPT for a certain mass range are 
the finite volume effects. There are two questions to consider: First of all, how does one know how 
large the zero mode contributions are, and how much do they affect the pion mass in the range 
of the fitting? Secondly, how large is the finite volume effect for the pion mass where the finite 
box is not larger than 4 times the Compton wavelength of the pion? We considered the zero mode 
contribution in Section ITTll where we examined the ratio of pion masses from the {P P) and {A4 A4) 
correlators for the 16'^ x 28 lattice and found that the ratio is basically unity within errors for the 
range of quark masses that we are concerned with. We now look at the same ratio in a smaller 
lattice 12^ x 28 {La = 2.4 fm in this case). They are plotted in Fig. Altogether with those from the 
16^ X 28 lattice. We see that the ratio for the 12^ x 28 lattice deviates from unity by two sigmas 
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for pion mass in the range of 0.5 to 0.8, or ~ 500 to 800 MeV. The central value of the ratio goes 

down to 98% at 771,^ ~ 300 MeV before it turns back to around unity for small pion mass due to 

the finite volume effect. This is a clear indication that the zero mode contribution to the {P P) 

propagator is visible even for pion mass as heavy as ~ 800 MeV. Most of the calculations compiled 

in Table mn have volumes similar to our 12^ x 28 lattice, therefore one needs to be concerned about 

the zero mode contamination which is a finite volume artifact. Since the zero mode contribution 
1 



varies as 



101,123, the effect of the zero mode contribution will be larger for calculations with 
still smaller lattice sizes. It would be useful to compare the pion masses from {P P) and (^4 A4) 
correlators to see how different they are in these cases. 
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FIG. 15: Ratio of pion masses as calculated from the {P P) and {A4 A4) correlators as a function of m,r 
from the {A4 A4) correlator for both the 16"^ x 28 and 12'^ x 28 lattices. 



To check the finite volume effect, we consider the pion mass from the {A4 A4) correlator which 
is less affacted by the zero mode contamination. We plot the ratio of pion masses between the 
12^ X 28 and 16'^ x 28 lattices in Fig. Elas a function of itIt^. We see that the ratio is basically unity 
for pion mass greater than ~ 300 MeV. Below this, the pion mass from the smaller volume, i.e. the 
12^ X 28 lattice, becomes higher. The worst deviation, at the smallest pion mass of 182(8) MeV, 
is 3.7%. For comparison, we also plot in Fig. ^Jthe corresponding ratio of from the {P P) 
correlator. Contrary to the {A4 A4) correlator, we see that the ratio has a dip between the m.^ 
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FIG. 16: Ratio of pion masses between the 12^ x 28 and 16"^ x 28 lattices for both the {A4 A4) and the 
{P P) correlators. 

range of 250 and 800 MeV, again reflecting the fact that the pion mass calculated from the {P P) 
correlator on the 12'^ x 28 lattice has observable zero mode contributions. 

To assess the finite volume correction, we first consider the leading finite volume correction 
calculated by replacing the infinite volume meson propagator with the finite volume counterpart. 
The correction [25| for is 

,2 



2tt 



(18) 



25l |. we can estimate the 



where /i^ = m^, - 771^/2 - mlo/2 = (871 MeV)^. With ami «; fi'^ 
percentage difference of between the 12^ x 28 and the 16'^ x 28 lattices for the smallest m,r 
182(8) MeV to be 

A{ml-^°°Pf{12) - A(mi-'°°P)2(16) 



0.13%. 



(19) 



This is much smaller than the Monte Carlo calculation of 7.4% derived from the (A4 A4) correlators 
in Fig. 1161 Next, we use the empirical dependence {41I to estimate the finite volume correction 
with 



mT,{L) = mT,{L = 00) + j^. 



(20) 



25 



For our lowest mT^{L = 16) = 182(8) MeV where the corresponding m-,^{L = 12) is 3.7% higher, 
we obtain the finite volume correction for the m.,^[L = 16) = 182(8) MeV to be 2.7% and the 
correction for m'^{L = 16) = (182(8)MeV)'^ to be 5.4%. Similarly, for the second lowest pion mass 
at 193(7) MeV, the correction is 3.6% for m^^L = 16). These are smaller than the statistical errors 
for the fitted low-energy parameters in Section HVI 

We are not in a position to fully address the continuum limit extrapolation, as we only have 



results at one lattice spacing. However, we know t 



lere is no 0(a) error due to chiral symmetry 
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and the O^a^) and O{'m?o?') errors are small [3, ^ We shall use the results from the small 
volume study in Ref. to gain some feeling for the magnitude of the O(a^) errors. We have 
calculated pion and rho masses on the 6^ x 12, 8^ x 16, and 10'^ x 20 lattices with the Wilson gauge 
action at /? = 5.7,5.85, and 6.0 respectively. The figures of mT^al \faa and mpa/^/aa (a is the 

n 

string tension) vs. cra^ in Ref. 8] show that both of them are quite flat at niT^/mp = 0.4, 0.5, and 
0.6. Since the lattice spacing at (3 = 5.7, as measured with the Sommer ro, is 0.171 fm which is 
very close to our present lattice whose lattice spacing as determined from tq is 0.175 fm, we shall 
use the /? = 5.7 data as the estimate for the present lattice. After extrapolating the pion mass 
from [3 = 5.7, 5.85, and 6.0 in Ref. 8j to the continuum limit, we find the central values of m,r at 
j3 = 5.7 are changed by 1 - 3 % for the cases with niT^/mp = 0.4,0.5, and 0.6. We thus tentatively 
suggest that our present pion masses are subject to the similar level of 0{o?) errors due to the 
continuum limit extrapolation. This, of course, needs to be checked with calculations at smaller 
lattice spacings than what we have now. 



VI. CHIRAL BEHAVIOR OF AND fp 



We show in Fig. El our data on renormalized /tt as obtained [lO| from 



A=Jim /"^^^^^^^~ • (21) 



We see that if one fits a linear curve in quark mass from 0.1085 (m^ = 0.4438(36) or 438(7) MeV) 
to 0.4 (ttItt = 0.8740(30) or 862(13) MeV) as is usually done in the literature, the fit gives f-j^ = 
0.0935(6) + 0.110(2) m with x^/NDF = 1.0. In this case, the data below m = 0.1085 shows a 



ow mass end. Even though in full 
in the form log(m^) with a 



hump with a one-sigma deviation from the linear fit at the 
QCD with dynamical fermions has a leading order chiral log 

known coefficient, it is absent in the quenched theory Q. We interpret the bump as a consequence 
of the quenched chiral log in the pion mass. Thus, we plot the the fir as a function of in Fig. 
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FIG. 17: /tt determined from the the pseudoscalar correlator Gpp{p — Q,t) in Eq. (|21|l as a function of the 
quark mass m. The dashed curve is a hnear fit to the high mass region, m = 0.1085 {ni^ — 0.4438(36) or 
438(7) MeV) to 0.4 {m„ = 0.8740(30) or 862(13) MeV). 



Eland fit with the following form 



fn = MO) + Ciml + C2mi. 



(22) 



As in the fit of m^, we fix the minimum pion mass at 0.1844(71) and fit the data to the above 
formula with variable ranges. The results of /7r(0) are plotted in Fig. ^1 as a function of the 
maximum pion mass w-^^^^.. We plot the same for Ci and C2/IO in Fig. [201 Even though there 
appear to be some changes in the region of m^^^ax between 0.3748(38) and 0.4840(32) for /7r(0), Ci 
and C2, the errors are such that one cannot be as certain as in the case of the chiral log 6. Thus, 
we shall quote the results for nij^^rnax = 0.3748(38) or 370(7) MeV. They are given in Table IVll 
The fit for in the range up to 

^TT^max — 0.3748(38) with the parameters in Table IVIl is plotted 
in Fig. ^1 as the solid lines for both the full figure and the insert. We also fit /^r with only the 
quadratic pion mass term for the large pion mass region from = 0.4438(36) or 438(7) MeV to 
0.8740(30) or 862(13) MeV) and obtain = 0.0952(6) + 0.0551(13)m2 . This fit is plotted as the 
dash line in Fig. 1181 

According to the quenched chiral perturbation theory ^2^. Ci = y^f^- From our data reported 
in Table rvn we obtain L5 = 3.6(1.5) x 10^'^. This is comparable to L5 = 2.5(5) x lO^'^ as obtained 



27 



0.16 I 

0.106 

0.15 - 

0.102 
^ ""^ " 0.098 
013 - 0.094 




0.09 - 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 

2 2 
a 

FIG. 18: as a function of m^. The dashed curve is a fit with /Tr(0) + Cim^ to the high mass region, 
= 0.4438(36) or 438(7) MeV to 0.8740(30) or 862(13) MeV). The sohd curve is the fit with Eq. ^ 
for the small mass region, m„ = 0.1844(71) or 182(8) MeV) to m^^^nax = 0.3748(38) or 370(7) MeV). The 
insert is a blowup for the small mass region. 

from the study of /tt with the clover fermion action on Wilson gauge configurations aX (3 = 5.7. 
But it is substantially larger than the phenomenological value of 1-4(5) x 10"'^ [4^. 

Using the experimental value fnifn-,^) = 92.4 MeV and m-,^ = 137 MeV to set the scale, we 
determine the lattice spacing a to be 0.200(3) fm. We see that fnim-jr) is higher than /7r(0) by 
-5%. 

The unrenormalized pseudoscalar decay constant of the pion is defined as 

/^ = (0|V^i75(l-I)/2)V'|vr(p = 0)) (23) 

TABLE VI: /^(O), Ci, and C2 as fitted from = 0.1844(71) or 182(8) MeV to = 0.3748(38) or 
370(7) MeV. 



MO) 




C2 


0.0905(23) 


0.150(63) 


-0.46(37) 3.6(1.5) x lO^^ 
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FIG. 19: /7t(0) as determined from the fit Eq. H22I) as a function of the maximum pion mass rn,r.,: 



From the chiral Ward identity, we have 



ZAdfj,Afj_ = 2ZmmZpP, 



(24) 



where A^^ = ^^7^75(1 — Z)/2)'0 and P = ^/5i75(l — Z?/2)i/^. Since = Zg^ and Zs = Zp due 
to the fact the scalar density S = -0(1 — D/2)ip and the pseudoscalar density P are in the same 
chiral multiplet, Zm and Zp cancel in Eq. H24|) . Therefore, from the relation of the on-shell matrix 
elements 



(0|ZA54^4|vr(^= 0)) = 2m(0|P|7r(p = 0)), 



(25) 



one obtains that 



u 



(26) 



/tt 2m' 

where f^^ is renormalized while fp and m are unrenormalized. The relation also holds when both 
fp and m are renormalized since Z^a — Zp^ . It is clear from Eq. (j26p that the chiral behavior of 
fp is the combination of those of m^/2m and /tt. 
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FIG. 20: Same as in Fig.Ulbut for Ci and C2/IO. 



VII. QUENCHED ttin 



We first plot the nucleoli mass as a function of in Fig. 1211 We see that it is fairly linear 
in all the way down to m,r ~ 0.4438(36) (438(7) MeV). Below that point, there appears to be 
a deviation from the linear behavior. Following the case of /tt, we shall first fit the mass only 
linear in ml in the range of ml = 0.196 - 0.763 and obtain mN = 0.998(21) + 0.932(74)m2 . This 
fit is represented by a dashed line in Fig. 1211 We see that the data (in the small m^^ region) in this 
case are systematically lower than the linear fit. 

According to the quenched xPT for the baryons Q|, the one-loop result predicts the presence 
of the non-analytic terms m-^^ and ml log(m7r), both of which are proportional to the quenched 5. 
Thus, we fit the nucleon data from the lowest pion mass, i.e. mT^^min = 182(8) MeV upward as done 
in previous sections with the form 



rriN = mAr(O) + Ci/2m^ + Ciml + Cuml log(^) + Cg/sm^, 



(27) 



We should mention again that the fit does not depend on the value of A^, since there is an analytic 
term in ml. We take A-,^ = 1.16 for this case. 

The fitted results of mAr(O) are plotted in Fig. [22l as a function of the maximum pion mass 



m 



■K,max 



of the fitted range. We see that it is rather stable in the range fitted. Similarly, the 
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FIG. 21: The nucleon mass as a function of mj. The dashed hne is a fit only hnear in for the range of 
from 438(7) MeV to 862(13) MeV. The solid hne is a fit with Eq. ^ for the range = 182(8) MeV 
to 316(6) MeV. 

TABLE VII: mAr(O), C1/2, Ci, Cil, and C3/2 as fitted from = 0.1844(71) (182(8) MeV) to 0.3200(43) 
(316(6) MeV). 



mAr(O) 


Cl/2 






C3/2 


0.957(43) 


-0.358(234) 


2.44(67) 


0.319(290) 


0.016(199) 



coefficients C1/2, Cil, Ci and C3/2 are plotted in Fig. We see that both C1/2 and Cil have 
rather smooth behavior in Tfiq^^jYidx from 0.3 to 0.55 and. approach, zero around, f^n^max — 

0.6 - 

0.7. On the other hand, Ci has more of a variation in the range mT^^max = 0.3 — 0.55. When 
nT'TTjinax reaches the region 0.6 — 0.7, the error bars do not touch those at fnir^max ~ 

0.3. We note 

that C3/2 is consistent with zero in the whole range of the fit. We have tried to include the next 
analytic term 46]. It turns out that mAr(0),Ci/2 5 and Ci are not changed outside their errors 
in this case. However, C3/2 and C2 with the term become ill-determined. They are consistent 
with zero but with large errors. Thus, the data seem to suggest that there is no evidence for a 
large term in the quenched nucleon mass. We quote the fit evaluated at rriT^^rnax = 0.3200(43) 
(316(6) MeV) in Table |VlII This fit is shown as the sohd curve in Fig. [21 for ml up to 0.4840(32) 
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FIG. 22: Fitted mAr(O) in Eq. H27|) as a function of rriT^^max of the fitting range. 

(ttItt = 478(8) MeV). Also shown in the figure as the dashed curve is the hnear fit from the higher 
pion mass range of = 0.4438(36) (438(7) MeV) to 0.8740(30) (862(13) MeV). 

We should mention that even though the enor . l„ge, . eons.teut with that predated 
by the one-loop xPT |45|, for which C1/2 = -^{D - 3F)'^6. Using 6 = 0.24(3) from our fit 
in Sec. IIVI it gives C1/2 = —0.37(5). This is consistent with the value in Table rvTll which is 
—0.358(234). Ci = 2.44(67) from the low mass fit is much larger than that of the high mass fit 
where Ci = 0.932(74). 

VIII. CONCLUSION 

The chiral extrapolation has always been a major challenge in lattice QCD calculations, pri- 
marily due to the fact that one does not know if the formula used to fit is correct in the range 
of the available data. The idea of using xPT to guide the fit is useful and appropriate except 
it is complicated by the fact that one does not know the range in where the one-loop formula is 
valid. In this work, we attempt to answer this question by fitting the data from our smallest pion 
mass (close to the physical value) upward and seeing where the low-energy parameters change. 
This is commensurate with the approach which expands around the small pion mass and 

momentum. The chiral log S from fitting the pion mass to the quenched one-loop chiral formula 
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FIG. 23: The same as in Fig. |55]for the coefficients Cil and C1/2 (upper figure) and for the coefficients Ci 
and C3/2 (lower figure). 



clearly demonstrates a sharp transition around 300 — 400 MeV in the pion mass. We take this 
as an evidence to reveal the upper limit for the applicability of the one-loop formula in quenched 
xPT. We have also fitted the pion mass with the cactus re-summed diagrams in quenched xPT 
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FIG. 24: The nucleoli mass as a function of for smaller pion masses. The solid line is a fit with Eq. (|27|l 
for the range = 0.1844(71) (182(8) MeV) to 0.3200(43) (316(6) MeV). The dashed line is a linear fit 
from a higher mass range of = 0.4438(36) (438(7) MeV) to 0.8740(30) (862(13) MeV). 

and the results suggest that its range of applicability extends beyond the one-loop formula. We 
should point out that fitting with the cactus re-summed diagrams is preferred over the one-loop 
formula since it gives a scale-independent S. However, not all the physical quantities have cactus 
re-summed results in the and in many cases, one would rely on one-loop formula to do the 
chiral fitting. Thus, it is essential to find out where the one-loop formula is valid. To be sure that 
this is not skewed by the finite volume effects, we compared our data on the 16^ x 28 lattice against 
those on the 12^ x 28 lattice and found that the finite volume errors are much smaller than the 
statistical errors. Although less certain due to larger errors, the analytic term in the mass and 
the non-leading chiral log in show variations in the pion mass range of 300 — 400 MeV, and the 
Ci coefficients change by a factor ~ 2 — 3 compared with those from fits involving exclusively pion 
masses beyond this range. Thus we tentatively conclude that the one-loop quenched xPT is valid 
for pion mass less than ~ 300 MeV. It is claimed that using a different regularization scheme such 

n 

as the finite range regularization 47] will extend the range of validity for the one-loop quenched 
xPT from that of the dimensional regularization, but it remains to be checked by lattice Monte 
Carlo data at low pion mass. 
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One of the valuable lessons learned in this study is that it seems futile to extrapolate to the chiral 
region from pion mass beyond ~ 300 MeV. We tried to fit the pion mass with the one-loop formula 
by extrapolating from ttItt = 862(13) MeV and found that, no matter how small pion masses one 
includes, 5 is inevitably different (by more than two sigmas) from that obtained with the data 
limited to <~ 250 MeV. One could try to model the functional form at higher masses, but it is 
not likely to succeed, because one does not know a priori how and when the modeled formula eases 
into that of the one-loop chiral xPT. Direct simulation at pion masses less than ~ 300 MeV is 
nearly impossible (or at the very least an extremely difficult task) for lattice calculations employing 
fermion actions which suffer from exceptional configurations and/or critical slowing down. An 
overlap-fermion calculation is numerically intensive, but with approximately fixed cost - its critical 
slowing down is rather mild The fact that it has chiral symmetry alleviates the worry that 
something may go wrong at small quark masses. Furthermore, the chiral behavior study is greatly 
enhanced by the incorporation of the multi-mass algorithm which facilitates the constrained curve- 
fitting method to discern the log structure. One does not know how different the chiral behavior 
in full QCD will be from that of the quenched approximation, but it is hard to imagine that there 
are different scales for the sea and valence quarks. If so, one will need to work harder to lower the 
dynamical quark masses such that the corresponding pion mass is below ~ 300 MeV. 
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